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\Q ■ We explore the application of approximation schemes from many body physics, 

f^^ ! including the Hartree-Fock method and random phase approximation (RPA), to 

T^ij- ' the problem of analyzing the low energy excitations of a polymer chain made up 

■^ ! of bosonic string bits. We accordingly obtain an expression for the rest tension 

'V ' To of the bosonic relativistic string in terms of the parameters characterizing the 

D , microscopic string bit dynamics. We first derive an exact connection between 

^ ■ 
L| ■ the string tension and a certain correlation function of the many-body string bit 

KA , system. This connection is made for an arbitrary interaction potential between 

c^ ' string bits and relies on an exact dipole sum rule. We then review an earlier 

calculation by Goldstone of the low energy excitations of a polymer chain using 

RPA. We assess the accuracy of the RPA by calculating the first order corrections. 

For this purpose we specialize to the unique scale invariant potential, namely an 

attractive delta function potential in two (transverse) dimensions. We find that 

the corrections are large, and discuss a method for summing the large terms. The 

corrections to this improved RPA are roughly 15%. 
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1. Introduction 

A decade and a half ago we proposed a microscopic theory of the relativistic 

[1-31 
string. At that time the main motivation was to seek a connection between 

string theory and QCD. From today's vantage point this motivation appears wrong- 
headed, but the model we proposed then still holds promise for a fundamental 
reformulation of string theory. Some of the same ideas have reemerged in recent 

years in the context of large A^ matrix models, which have been the subject of 

[4—81 [91 

vigorous investigation. In 1991 we sketched why the wee parton model of 
Ref . [3] is an attractive starting point for a fundamental formulation of string theory. 

In this article we begin a serious study of the detailed relationship of the 
microscopic physics of such a string bit model to the low energy properties of 
string. Presumably, the microscopic details of the model would only be apparent 
at the Planck scale. For the purposes of our present analysis we imagine that the 
string coupling is small, i.e. the string mass scale ^/Jq « Mpianck- Even though 
this is probably an unrealistic assumption, it is useful because it is a limit which 
puts highly nontrivial constraints on the string bit model. In particular, it is a limit 
in which the graviton, a composite of string bits in our model, must satisfy all of 
the low energy theorems enjoyed by Einstein's general relativity. These constraints 
are guaranteed, of course, once it is demonstrated that the low energy properties 
of our string bit system are precisely those of the relativistic string. 

The work of Ref. [3] established the microscopic/macroscopic connection be- 
tween string bits and the bosonic string in the context of the Random Phase Ap- 
proximation (RPA). It can be argued that the corrections to RPA can be lumped 
into a renormalization of the string tension, and so are not important as far as 
the stringy properties of the low energy theory are concerned. However the actual 
size of the corrections to RPA was not determined in Ref. [3]. Since the actual 
microscopic details of the string bits may be of eventual interest {e.g. in resolving 
the short distance properties of gravity), we believe that it is important to assess 
the validity of RPA quantitatively, and this is our goal here. Moreover, the de- 



tails of the micro/macro connection are probably essential for a proper string bit 
understanding of the superstring. For example, must we require an exact super- 
symmetry at the level of string bits, or is supersymmetry only present in the low 
energy theory? We don't answer that question here, but it is clearly an important 
one, and the answer is surely going to require a more extensive treatment than 
that given in Ref. [3]. 

In the following sections we provide this more detailed study of bosonic string 
bits. In Section 2, we present the dynamics of a polymer of string bits in Gold- 
stone's formulation of the many body problem. This formulation provides an 
elegant presentation of various "improved" versions of perturbation theory, such 
as Hartree-Fock methods and RPA, in terms of Feynman diagrams. In Section 3, 
we use this diagrammatic method to define certain "irreducible" correlation func- 
tions, for which exact low frequency theorems can be derived. These theorems 
are then used to show the universal presence of harmonic low energy excitations, 
which reproduce the usual spectrum of the relativistic string. Thus the string 
tension will be related to the value a certain zero-frequency correlation function 
J dt {G\ Tx'^{t)x''^{0) \G). This holds even for string bit interaction potentials that 
are not differentiable or even finite at minimum separation. A case in point is the 
potential V(x) = — Ao5(x), which has the nice property of being scale invariant 
in two dimensions — precisely the case for the transverse space of four-dimensional 
space-time. But evaluation of the relevant correlation function requires an exact 
solution of the string bit system, which we don't have except in the case of a har- 
monic oscillator string bit potential. Our string bit model involves a short range, 
scale invariant potential. This is where we must resort to approximations. In Sec- 
tion 4, we review Goldstone's application of RPA to the string bit polymer problem 
and then proceed to calculate the first order corrections to RPA. We find that the 
diagrams which have the character of "self-energy" corrections are in fact quite 
large. Fortunately the structure of these large corrections is simple enough that 
they can be summed to all orders, which we do in Section 5. The upshot is that 
the large corrections have the effect of renormalizing the ionization energy of the 



polymer, and once this is taken into account the remaining corrections are small, 
of order 15%. Finally we gather some conclusions in Section 6. 

2. Chain Dynamics in Formalism of Second Quantization 

The model of string bits proposed in Ref.[3,9] leads to the relativistic string 
in light-cone gauge. Although the string formed from string bits moves in D 
space-time dimensions (described by light-cone coordinates 






where d = D—lis the spatial dimension), the string bits move only in the transverse 
space, with coordinates x: in fact they enjoy a Newtonian {i.e. Galilean invariant) 
dynamics in this transverse space. The x^ = t coordinate is a Newtonian time 
for the string bits, but the longitudinal dimension, corresponding to x^ and its 
conjugate momentum P~^, is nonexistent for them. It only emerges with string 
formation: bit number is a conserved quantity in string bit dynamics, and for a 
string with a large number of bits it becomes an effectively continuous variable, 
which can be consistently interpreted as a variable P+. In this interpretation, we 
can think of each string bit as carrying a tiny fixed amount of P+ = e. Then 
a string containing M string bits would carry P+ = Me. The fact that the P~ 
of the string should vary as 1/P^, is then understood as the 1/M dependence 
of low energy excitations of large systems. The hamiltonian describing string bit 
dynamics will accordingly be identified as the operator eP^ for strings with an 
infinite number of string bits. 



The string bit dynamics proposed in Ref.[3,9] exploits 't Hooft's ideas on the 

[121 

1/A^ expansion. We postulated that the string bit annihilation operator is an 
N X N matrix variable Oj^ (x) and denoted the string bit creation operator by 



a/(x) = a/(x)t. Then 

In the large N hmit, one can show that the singlet operators 



N 



A(xi, ■ ■ ■ , xm) = ( — J tr{a(xi) ■ ■ ■ a{xM)} 



[131 

behave as creation operators for closed polymers of string bits. Note here that 
Bose statistics for the string bits implies cyclic symmetry for the A's. This sym- 
metry accounts for the well-known Lq = Lq constraint in the eventual low energy 
string dynamics. 

A candidate dynamics for bosonic string bits is given by the Hamiltonian 

H =eP' = / rfx- tr Va(x) • Va(x) 

^ 2iV / ^^^y^^^ - y) tr : (a(x)a(x) - a(x)a(x))(a(y)a(y) - a(y)a(y)) : 

(2.1) 
where V should be positive (repulsive) for stability. Evaluating the action of H on 
a state 

1 \ M/2 



leads to 



, tr{a(xi) ■ ■ ■ a(xM)}^M(xi, ■ ■ ■ , xm) |0) 



H m = I Ef (-ViV2 - V(x,+i - x,))^) 
1 ^ - 

-|- Other Terms of Order (1/N). 

(2.2) 

The displayed 1/A^ term describes the dissociation of a closed polymer into two 

closed polymers, so 1/A^ can be identified with the string coupling constant. In 



the limit that it vanishes (A^ -^ oo), the above equation describes a stable polymer 
provided that the nearest neighbor potential —V is sufficiently attractive to bind 
a pair of string bits. We shall review the argument below. 

Our aim in this article is to go as far as we can toward the solution of this large 
A^ limit of the string bit system. Referring to (2.2) we see that the energy eigen- 
states in this limit are determined by the solution of the following nonrelativistic 
many body Schrodinger equation: 

M 

^(-Vi2/2 - V(x,+i - x,))^(xi, . . . , xm) = ^t^(xi, . . . , xm). (2.3) 

For purposes of keeping this article reasonably self-contained, the rest of this section 
will be devoted to a recapitulation of the main results already reported in Ref.[3]. 
New developments will be reserved for the subsequent sections. 

It is first convenient to use Galilean invariance to separate the motion of the 
center of mass. So define 

^V (2.4) 

where xq = xjvf and XAf_|_i = xi. Note that although we have defined a y,. for each 
r, only M — 1 of them are independent: they satisfy the constraint 

1 ^ 
M ^^"^ 

r—l 

Choosing y2, . . . ,yM as independent variables, and denoting (l/i)(9/9X = P and 
{l/i)d/dyr = Pr, the many body Hamiltonian corresponding to (2.3) becomes 

p2 M M-l M 

^ ^OM + 5Z P^^ ^ 5Z P^ ■ P^+1 ~ Yl ^(y^) ~ V(-y2 jm) 

^^^^ r=2 r=2 r=2 (2.5) 

p2 
= — + h'. 

2M 
The first term on the r.h.s., giving the energy of center of mass motion, is precisely 



what is needed for the consistent interpretation of Me as P+ and h/e as P . The 
remainder of the r.h.s. h' should have the interpretation (for M -^ oo) as ]v(^ J2M 
where M. is the rest mass operator for the string. Since there will inevitably be a 
bulk contribution to the eigenvalues of h! proportional to the number of string bits 
M, we must subtract this contribution before the \jM behavior is revealed. This 
subtraction is easily accomplished by including a counterterm in the original string 
bit Hamiltonian of the form /iq / rfxtr a(x)a(x). Note that in the low energy theory 
such a term is just a multiple of P+. The remainder of this paper is devoted entirely 
to gaining an understanding of the eigenspectrum of h! especially the excitation 
energies of order (1/M) above the ground state energy. 

We first use variational arguments to establish a simple criterion for the stabil- 
ity of the string bit polymers. Here we are assessing stability against the breaking of 
nearest neighbor bonds, i.e. the dissociation into a number of open polymers with 
all unbroken bonds retaining their arrangement. Dissociation of a closed polymer 
into two smaller closed polymers by bond rearrangement is another matter inti- 
mately associated with the presence or absence of tachyons in the low energy string 
theory. Such decays are vetoed by the N ^ oo limit corresponding to zero string 
coupling constant. 

We assume that the potential —V is negative and vanishes at large distances. 
The zero of energy corresponds to all the string bits infinitely dispersed from one 
another. Notice that deleting the last term in (2.5) leaves the Hamiltonian for an 
open polymer with the same number of bits. Since the deleted term is negative, 
this shows that the closed polymer ground state energy is strictly less than the open 
polymer ground state energy. Thus stability of the open polymer implies that of 
the closed polymer. Next consider the dissociation of an open polymer with M 
bits into two open polymers with Mi and M — M\ bits respectively. Because of 
the nearest neighbor interaction structure, the internal Hamiltonian for an open 
polymer can be broken down as: 

^Open,M = ^Open,Mi+^Open,M-Mi+PMi + l~"'^(yMi+l)-PMrPMi+l-PMi+rPMi+2, 



where the variables labelled by Mi + 1 do not appear in either of the open polymer 
Hamiltonians on the r.h.s. Thus applying the variational principle with a trial 
wave function given by ^Trial = ^Mi^M-Mi^2^(yMi+i), where ^f is the exact 
ground state wave function for a polymer with K bits, gives 



TpG ^ TpG I TpG I TpG 

-^Open,M ^ -^Open,Mi + -^Open,M-Mi + ^2 ■ 



Here one uses the fact that (pr) = in any of the ground state wave functions in the 
trial. If the nearest neighbor potential binds two bits, E2 < and the ionization 
energy required for dissociation into two open polymers is at least B = — -^2 • 

Now for large M, the ground state energy of an open polymer should have the 
behavior 

Inserting this form into the variational estimate gives the rigorous bound f3 > 
—E2 = B. More precisely we see that (3 has the interpretation as the ionization 
energy required to break a long open polymer into two smaller but still long open 
polymers. This ionization energy is just the energy cost of breaking a bond. Since 
the transformation of a closed polymer to an open polymer requires breaking pre- 
cisely one bond, and for very long polymers the energy cost of this break can only 
differ from that of breaking an open polymer into two open polymers by energies 
of order 1/M, we can conclude that the ground state energy of a closed polymer 
should have the behavior 

^Closed,M~«^ + 77M. 

The coefficients of the 1/M terms are, of course, different in the open and closed 
cases. It is crucial for the interpretation of h as eP~ that there is no constant (3 
term for the closed polymer ground state energy, since only the bulk a term can 
be removed by a counterterm in the string bit Hamiltonian. 



The important parameter that characterizes the stringy properties of long poly- 
mers is the string rest tension Tq. It can be inferred from the low lying level spacing 
of the polymers. Eigenvalues of h' of order 1/M above the lowest eigenvalue will 
correspond to values of P^ behaving as 1/P+. The coefficient of 1/M is thus 
2 X A^^. To calculate Tq it is simplest and sufficient to study long open polymers, 
and we shall make this restriction in what follows. 

We only remark that the best way to apply the methods to the closed case is 
to consider a related system to (2.5), in which one new degree of freedom (yi, pi) 
is introduced, and the Hamiltonian is taken to be 

p2 M M M 

r=l r=l r=l 

To apply solutions of this system to (2.5), notice that the quantity Y, defined 
below (2.4), commutes with h. One can therefore look for eigenfunctions of h of 
the form 5(Y)-?/'(y2, . . . , Ym); whence ip must be an eigenfunction of (2.5). Thus the 
eigenstates of (2.5) are the subset of eigenstates of h which satisfy the constraint 
Y = 0. 

Henceforth we restrict attention to the open polymer hamiltonian with center 
of mass motion removed 

M M-1 M 

r=2 r=2 r=2 

Our variational arguments have already established that the ground state of this 
system is a discrete bound state lying at least an amount B below the continuum. 
We must next look for discrete levels lying 0(1/M) above the ground level. For 
this purpose, we apply many body methods developed by Goldstone. ' 

Goldstone employs a second quantization formalism in which a fermionic field 
ipriy) is introduced for each bond r = 2, . . . , M in the open polymer chain. Each 
bond number operator Mj. = J d'yijjl{y)tjjr{'y) will be conserved by the second 



quantized dynamics and so can be consistently set to unity. In this sector the 
second quantized Hamihonian 



r=2 
M-1 

-J2 J dy4{y)H'^)My) J dz-4^,{z){-tV)A+i{^ 



r=2 

is completely equivalent to (2.6). This Hamiltonian can now be analyzed in the 
interaction picture using Feynman diagrams in which the propagator is determined 
by the quadratic term and the vertex is determined by the quartic term. The con- 
ceptual advantage of this formalism is that the zeroth order ground state, namely 
the state in which each bond is in the ground state of its one body potential, is 
simply expressed as a boundary condition on the propagator, just as in the Dirac 
hole theory of positrons. 

Expanding ■?/',■ in eigenf unctions of the one body hamiltonian p^ — V(y) 



^r(y) = ^ &m0m(y) + &S0o(y), 



m>0 

where m = Q labels the ground state and m > labels the excited states, the 
zeroth order ground state is specified by 

hi, |0) = for m > feg^ |0) = 0. 

Then the interaction picture propagator, represented by the first diagram in Fig. 1, 
is given by 

=5rs{e{t) Y, 0m(y)0;.(z)e-^*^™ - ^(-t)0o(y)0a(z)e-'*^n 

m>0 (2.8) 

d^ -tcut i \^ -'^MyWmi^) , -i0o(y)0o(z) 






27r -^^ Em — uj — ie En — uj + ie 

' m>0 
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Fig. 1. Feyninan rules for a string bit chain. 
The vertex function represented by the second diagram in Fig. 1 is given by 

i6r^sA2S2{^r^,r2 + l + (^ri.ra-l) (-^Vl)(-iV2) (2.9) 

where the gradient operators act on the first arguments of the respective propaga- 
tors directed into the vertex, and a Kronecker delta with either index equal to 1 
or M + 1 is understood to vanish. 

In order to apply time dependent perturbation theory to the calculation of 
ground state expectation values of time ordered products of operators like y^ or p^, 
we simply use Dyson's formula after translating to the second quantized formalism. 
Thus we have 

/I d 

so, for example, 

{G\T[y'At)yTm\G) = 

(0| T[f dyjj'rHy, t)y>'ijUy, t) I dy^plHz, 0)^-^,^(z, 0)e-'^-^ '^'^^^^'^] |0) 

(2.10) 
where Hj{t') is just the quartic term of (2.7) transformed to interaction picture. 
As usual the denominator in (2.10) simply cancels all disconnected bubble dia- 

11 



grams and correctly normalizes the matrix element after the infinite time integral 
projects onto the true ground state of the many body system. Information on the 
spectrum of the exact Hamiltonian can be inferred from the time dependence of 
such correlation functions, or equivalently from the pole locations in their time 
Fourier transforms. 

By applying the Random Phase Approximation corresponding to summing the 
bubble diagrams shown in Fig. 2, Goldstone demonstrated the existence of order 
1/M excitations in this approximation and gave an approximate formula for the 
string tension. In the following section we put these insights into an all orders 
context, and so obtain an exact formula for the string tension. 



Fig. 2. Diagrams summed in the random phase approximation 



3. Low Energy Excitations: The String Tension 
Consider the system 



M M-1 

^Open = X1[P^ + ^(yr)] " J^ Pr ' Pr+1 
r=2 r=2 

Our aim is to calculate the low lying excitations of such a system under the as- 
sumption that the system binds together, i.e. that there is a discrete ground state 
with a mass gap. We shall examine the poles in the time Fourier transforms of 

12 



'tut 



iuot 



various correlation functions: 

{G\T[yi{t)yim\G) ^ - j ^S%{u^)e 
{G\T[pl{t)pim\G) ^- J ^^g%{uj)e 
{G\ T[yi.{t)p^M] \G)^- j ^^.r^(^)e-^^^ 

{G\ nmyim i^) = - / ^^^rV^i^y-''^'- 

Using time translation invariance, we can transfer the time argument to the other 
member of the time ordered product, which leads to the identities 

Srs{-(^) = Ssr{i0) 

Using the Heisenberg equations of motion we easily derive the "Ward identities" : 

I {G\T[yl{t)yim \G) = {G\T[{2p^,{t) ' I^r^i^t) - I^r-iUMm \G) 

I {G\T[yl(t)TPM] \G) = ^S^j6rsS{t) + (G| T[(2p^(t) - p^+i(t) -p^_i(t))p,(0)] \G) 

Taking Fourier transforms of the Ward identities, defining k'^g = 25rs—Srs+i—Srs~i, 
leads to 

-liuS'^iiu) = k^r^y{uj) (3.1) 

-lur^iu) = 6yl + eg'^{u) (3.2) 

iurjy{uj)=5i,l + g'^{oo)e (3.3) 

luS'^uj) = rjp{uj)k^. (3.4) 

The bond labels r, s have been suppressed in writing these equations, i. e. each of 
the quantities S, J-'py, and Q is an (M — 1) x (M — 1) matrix, and the multiplications 
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indicated in these equations are matrix multiplications. Corresponding to £, J-'py, 
and Q, we shall introduce "irreducible parts" Hyy, ^py, ^yp, ^pp- Referring to the 
Feynman diagram expansion of Section 2, these quantities are defined as the sums 
of all diagrams, contributing to the respective correlation function, which cannot 
be cut in two by cutting a single dashed line. See Fig. 1. Then we can write an 
equation for each of the quantities S, J-'py, and Q as follows: 

^^'M = n;-^M[<5,fc + (A;2-2)^^-^] (3.5) 

J^iliu;) = Ul{iuj)[6,k + ik' - 2)g^'] (3.6) 

•^SM = [5.. + ^^^'(A:'-2)]ng (3.7) 

S'Huj) = U% + J^g{k^-2)U^. (3.8) 

This set of equations gives a purely algebraic determination of the irreducible parts 
in terms of the full correlation functions. 

Next we wish to use the Ward Identities to derive relationships among the 
various II's. First use (3.2) and (3.6) to eliminate J-'yp: 

5^JI + k^g'\uj) = -tuY['y\{u)[5kj + (A;2 - 2)g% (3.9) 

and then (3.5) can be used to rewrite the l.h.s. of (3.9): 

5^jI + k^g'\uj) = 6^jl + (A;2 - 2)g'^ + 2Uil[6kj + {k^ - 2)g^^] 
= [5,k + 2W^^]{5k, + {k^-2)g^^). 

Putting this into (3.9) gives, 

[5,k + 2n;^ + iul]^y%]{5kj + (A;2 - 2)^^-') = 0. 
Now provided 5kj + {k'^ — 2)g -^ is invertible, we obtain our first relation: 

-tuUi{ = 6ij+2Ui{. (3.10) 

Next, the identities g{—uj) = g'^{uj) and J^py{uj) = J^^ J—uj) yield the transposed 
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relationship: 

tujU'^y{uj) = 6ijI + 2U%{u). (3.11) 

The last relation we seek is between Uyy and Upp. First insert (3.4) in (3.8) to 
obtain 

r^(c.) = u% + ^c.r^(c.)^ng(c.). 

Then use (3.11) to elimintate Hpy, 

£^\uj) = 11% + £^^{u;)^^i5,kl + 2ng(c.)), 

Leading to 

U%{uj) = rH^)[pjkI - ^2ng(c.)]. (3.12) 

(3.2) and (3.4) together give 



uj^s'^l^ = Siji + k^g'^{uj) 



which, inserted into (3.12), yields 

u^ufyiu) = 2[6ijl + k^g'^{u)][6jkl - {k" - 2)n^J(c.)] 

= 25ikl + 2k'g'\uj) - 2(A;2 - 2)n;^(cc;) - 2k'g'^ {uj){k' - 2)^. 

(3.13) 

Now putting {uo -^ —uo) and taking the transpose of (3.5) shows that 

^^^-n;^ = ^^^'(A;2-2)ng. 

Using this result in the last term of (3.13) gives, after some simplification, the last 

15 



desired relationship: 



or, rearranging a bit, 



;24^,M = 25,fc/ + 4n;^( 



n;» = -^5,fcj + ^nj^^^H. (3.14) 

The imphcations of the Ward identities on the irreducible parts are summarized in 
(3.10), (3.11), (3.14). 

Armed with these results, we now turn to their implications for the low lying 
excitations of the polymer system. This information is given by the location of poles 
in Q, or by the zero eigenvalues of its inverse, which, for cj ^ 0, is proportional to 

Sjk - Ilf,{uj){k^ - 2) = 6,k + Uk' - 2)6,k - ^ng(A;2 - 2) 

2 (3.15) 

= -A;2<5,fc-^ng(c.)(A;2-2). 

The matrix k'^ has eigenvalues of 0(1/M^) as M ^ oo, so there are eigenfrequen- 
cies u'^ of the same order. For these we may neglect /c^ in the second term of (3.15) 
and evaluate Uyy at cu = 0. We therefore must search for zero eigenvalues of the 
matrix 

u'^UfyiO) + kHik (3.16) 

as M ^ oo. It is interesting and useful to note that by virtue of the Ward 
identities, £ (0) = 11*^(0). Thus there is no a priori requirement to separate 
out the irreducible part: at zero frequency the reducible parts vanish identically. 
Finally, by rotational invariance, we can write 



SO that (3.17) can be simplified to 

co^UyyiO) + k^ (3.17) 

16 



To examine the large M dependence of (3.15) go to normal modes: 



2 v;^ mnr nns 

i^mn = 7 Sm Urs Sm 

r,s 

'^mn = '^mn2(l — COS — ), 



after which the r.h.s. of (3.15) becomes 



djkOmn[^ - COS — ) + —Iimn[^)Ojk COS — . 

To get the spectrum of the relativistic string as M -^ oo, we should find that 
uj oc n/M, which requires that nmn(O) must be diagonal and independent of m for 
small m/M or n/M . The condition for this to be true is that, away from r = 2, M, 
rij-s — > sufficiently rapidly as | r — s |^ oo. 

To see this, calculate 



2 v;^^ . rrniir + s) . mis 
sm sm 



^rnn " ,^ /^ ^r^ 



M ^^ ^' M M 

r,s 



For r, s away from the endpoints, Urs should be only a function of r — s, so denoting 

it as n(r — s), we have 



2 v^T-T/ nns 17171 s mirr rmrs mirr, 
llrnn = — > Il(r) sm (sm cos h cos sm 

s,r 



/■^-^ mnr \ fsr^-n/ mvrrx 2 ^-^ nirs mirs 

Omni/ cos— — — n(r) + > n(r ) sm — — — -— > sm — -— cos 

r 
<5mnEn(r)] 



M 

r s 



as M -^ OO, provided W (r) vanishes rapidly enough as r -^ oo. This is cer- 
tainly true to any finite order of perturbation theory. Making this assumption 
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nonperturbatively, we have, 



mir, , . o mir m it 



a;2 Y^ n(r) ^ -2(1 - cos — -) = -4 sin^ 



M^ 2M (M)2' 

Since to is to be identified with eP~, and P^ = eM, the (mass) = 2P^P^ is 



given by 217171/ y^— ^^. n(r), by definition, the coefficient of m in this formula is 
l/a' = 2txTq: 



To = , ^ ^ a' = ^ /-yn(r) (3.18) 

V/-E.n(r) 2vrPW Z. ^ ^ ^ ^ 



4. Random Phase Approximation and First Corrections 
The Random Phase Approximation (RPA) is to calculate n^i replacing H by 

ngpA,..(o) = 5..ff^'(o) 

where 

CO 

tf^'(O) = -i j dt{0 I T[y\t)y^{0)] \ 0) 
—00 

corresponding to the one body problem 

^Ibody = P^ - V(y) 

Putting in a complete set of states: 
00 
tf^'(O) = -^ / dt[e{t) J2 (0 I ^'(0) \n){n\ y^{0) |0) e-*(^"-^")*+ 



-2E 



t) 5^(0 I l/^'(0) I n)(n I y\0) \ o)e+*(^"-^°)* 
(0|y|n)-(n|y |0) f 6,/ 
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1/2 
RPA J V^ yOn • ynO I ^T 



From which, the Regge slope is 



c^RPA 




1/2 



YnO 






fl4l 

a resuh first obtained by Goldstone. 

On the face of it, it would seem that the validity of RPA depends on the 
smallness of the perturbation. That this is not necessarily so is shown by the 
fact that for —V = Ky^ the RPA gives the exact result for the string tension, 
even though the ground state energy receives finite corrections. This is by no 
means transparent in terms of diagrams: the corrections to the string tension 
cancel, apparently miraculously. One can expect that potentials that are close to 
quadratic near the minimum will lead to a quite accurate RPA. 

However, probably the most interesting candidate string bit potential is a local 
one, V(y) = Ao5(y). In two transverse dimensions (corresponding to four dimen- 
sional space-time!), this is a naively scale invariant potential, although ultraviolet 
divergences give logarithmic scaling violations. Scale invariance is a desirable prop- 

[151 

erty of the microscopic potential since, as Gross and Mende and others have 
emphasized, it should be enjoyed by string theory at short distances. Even apart 
from its appeal for string bit dynamics, this potential is drastically different from a 
quadratic one, which makes it a good test case for assessing the accuracy of RPA. 
Thus in this section we consider corrections to RPA for this special potential. 

We begin by reviewing the basic features of the one body problem represented 

by 

/iibody = P^ - Ao5(y), (4.1) 

as it determines the zeroth order in time dependent perturbation theory. First 
consider the bound state problem. For an energy eigenvalue E = —B, B > 0, the 
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momentum space wave function is easily shown to be given by 

Ao^(O) 

Since ^{0) = J (ip0(p)/(27r)^, the quantization condition for B is obtained by 
integrating both sides over p: 

_ . f dp 1 

In two dimensions the integral diverges in the ultraviolet, which requires asymptotic 
freedom (Aq -^ 0) in order for B to be finite. With an ultraviolet cutoff A on the 
momentum integral, we find 

B = A^e'^^/^\ (4.2) 

[31 
that is we have dimensional transmutation, in which the dimensionless Aq is 

traded for the dimensionful B. In the limit A ^ oo, '?/'(0) is infinite but Ao'?/'(0) is 

finite. The normalized wave function is 

^(y) = ^i3(|y|) ' ' 



TT y 27r p^ + 5 .^ g. 

-ifodylv^), 

TT 

where Kq is the zero order modified Bessel function of the third kind. 

Clearly this zero angular momentum state is the only discrete state. The rest 
of the s-wave spectrum is positive energy continuum. For each energy E = /c^, 
there is one state whose wave function is that unique linear combination of Bessel 
functions Jo(/i;|y|) and Fo(^|y|) = ^o(^|y|) which is orthogonal to (4.3): 



^fc(|y|) 



,^°(^I^I) + MW^^°^^I^I) 



(4.4) 



where we have normalized R^ according to 

/27r 
dyRk{\y\)Rk'{\y\) = ^S{k-k'). (4.5) 

The radial wave functions of the states with nonzero angular momentum Lz = m 
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vanish at the origin and so coincide with the corresponding free wave functions, 
namely i?-(|y|) = J^(A;|y|). 

From the above information we can easily construct the Green function 



GEiy, z) = (y| T T^— |z) 

/^ibody - rj - %e 



The s-wave contribution to Ge is just 



i?B(|y|)i?B(|z|) , rkdkRki\y\)Rki\z\) 





^^ ^^'^^" -B-E-ie +/ 27r k^ - E - le 



The non s-wave contribution is just the free Green function minus its s-wave. Thus 
we have 

/• rfp e^P-(y-^) i?g(|y|)i?g(|z|) 

^ °^kdkRk{\y\)Rk{\^\)-Uk\y\)Jm^\) ^^'^^' 



2n k^ -E-ie 



By working directly with the two dimensional partial differential equation that 
determines Ge, it is straightforward to show that an alternate representation is 



GE{y,z 



(27r)2 p^-E-te ^^^^^ 

Apart from the distribution of ie's, Geo is just the Fourier transform of the propa- 
gator (2.8) needed for developing time dependent perturbation theory. To get the 
right boundary conditions for the propagator, one simply has to change the —ie in 
the middle term of (4.6) to a +ie. To keep track of this change, it is convenient to 
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define a reduced Green function Ge wliicli doesn't liave tlie bound state pole pole: 



Ge = Ge 



RB{\y\)RB{\A) 
-B-E-ie ' 



Then the propagator (2.8) can be expressed as 



T[^Uy,tW{^,o)]\o) = -tSrs I ^e 



= — i6r 



du 



'iiot 



Gu 



duo 

2^' 



e-*^^a 



RB{\y\)RB{\A) 
—B — uo + ie 



(4. 



With this explicit representation for the propagator we are now in a position to 
write down and evaluate the expressions assigned to the first few diagrams in time 
dependent perturbation theory. 

Focusing on the string tension, we have seen in the previous section that we 
need to compute Tlyy at zero frequency. To lowest order in perturbation theory, we 
must evaluate only the first diagram of Fig. 2. (The remaining diagrams in Fig. 2 
are reducible, and in any case vanish at zero frequency.) In next order we have the 
diagrams of Fig. 3. To this order irreducible diagrams with the external sources 
attached to different lines vanish by rotational invariance. 



Fig. 3. First order corrections to RPA. 
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Turning first to the lowest order diagram, the Feynman rules assign to this 
diagram the expression 

\ —B — oo + le J \ —B — uo + le J 

The prefactor phases have the following origins: the first —i from the definition of 
n, the {—tf' from the two propagators, and the (—1) from the closed fermion loop. 
The uj integral can immediately be done by closing the contour in the upper half 
uo plane, picking up the pole at u; = —B + ie. (The singularities from G are all in 
the lower half plane.) The result is 



-25rs j dydzRBiy)y'G^Biy, z)/i?B(z), 



which is identical to the result of the previous section. Because Rb is an s-wave 
state, the presence of y^ and z cause all s-wave contributions to G to decouple, 
so the latter can be replaced with the free Green function. Working in momentum 
space we get 



oo 



B 

This approximation for Ilyy is, as we have discussed, the random phase approxi- 
mation, giving a string tension of Tq = \I2txol ~ Tq = i?v3- 

Corrections to RPA are given by the three diagrams in Fig. 3. The first diagram 
represents the expression 

dijj\ duJ2 du^ 



25r-ii(—i) 



-1) i / dydzdudvdxdMV I 



27r 27r 27r 
y^G^, (y , u) (-2 V;r)G'a;i (u, z)/^^, (z, w){-iV'i)G^, (v, y) 

(-iV")G',;,_c^2+^3(x, w)(-iV;^)G'a;3(w, x) 

The factor of two comes from the two possible terms in the vertex function (2.9); 
this assumes that r, s are well away from 2 and M . To begin reducing this ex- 
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pression first note that the uj^ integral can be closed in the upper half plane where 
there is only a simple pole at 0^3 = —B + it. Thus 



27r 



cixdw(-iV")G'cji_,^2+cj3(x, w)(-iV^)G'(^3(w, x) 



i / rfxrfw[i?B(x)(-iVS)G'a;i-c^2-B(x,w)(-iV;;;)i?B(w) + 



+ 



i?B(w)(-2V:^)G^,„^,_B(w,x)(-2VS)i?B(x)] 

1 






-iB6r. 



dt- 



t-B 



2 J ~ t2 

B 



p'^ + B — uji + UJ2 — 2ie p^ + B — UJ2 + uji ~ 2ie 
1 1 



+ 



t — uji + UJ2 — ^.ie t — UJ2 + uJi — 2ie 



(4.9) 



This same integral appears in the other two diagrams, so we shall use this result 
again. 

The first diagram has now been reduced to 



2i5rs / dydzdudv 



duo I duj2 
27r 27r 



(4.10) 



-iB6r^ 



dt 



t-B 



+ 



t — uji + UJ2 — 2ie t — UJ2 + uji — 2ie 



B 



Now consider separately the two terms in the last factor. In the first term close 
the uji integral in the upper half plane, but in the second close the 002 integral in 
the upper half plane. In each case there is only a simple pole at ui = —B + ie, 
respectively 002 = —B + ie. The potential double poles at these locations don't 
contribute because they correspond to two s-wave states separated by a vector 
operator. (This will not apply to the other two diagrams!) After this the remaining 
CO integral is done by again closing in the upper half plane. Altogether there are 
four separate contributions, each of which can be thought of as a ground state 
expectation value in the one body problem. The results of the steps described so 
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far can be represented, in a self-explanatory notation, as 



t-B. 



-B6rs / dt- 2 
B 



{B\y 



p^ + ±5 p^ + i3 + t 



p- 15) 



- 7 (^1 y'^^p"" \B) {B\ y'^^f^ \B) 

t Y* + B p^ + B 



{B\f 



^ y'G^t-BP'^Tr^y' \B) 



.k 



p2 + 5 + t 



p^ + B' 



V^\^-^y'\^)^^\^-^ 



y'\B) 



+ {B\p'^ 



^ y'G-t-BP'^-^y'lB) 



p^ + B + t 



p' + B' 



1 



\(^\p-^y'\^)(^\p-^.^B 



y' \B) 



+ {B\y' 



p2 + 5 



P (^-t-By 



1 



P 



2 + fi + t 



p^ \B)] 



UB\y'^^P'^\B){B\y'^^ 
t p^ + i5 p^ + B 



P^ \B)] 



(4.11) 



where we have replaced the first and last Green functions in each term by the 
free ones, since the s-wave parts decouple in these locations. We postpone explicit 
evaluation until we have reduced the remaining diagrams. In addition, e may be 
safely set to zero in all of the terms in (4.11). 

We present the expressions for the second and third diagrams of Fig. 3 after 
the first simplification step analogous to (4.10). The second diagram leads to 



— 2i5rs I dydzdudv 



duji duj2 
l/JG^,(y,z)^'=G,,(z,u)(-^v:r)G^,(u,v)(-^v;r)G^i(v,y) 

oo 

-iB [ t- B \ 1 1 

at- 



(4.12) 



t2 



B 



t — cji + c<J2 — 2?e t — CJ2 + ^1 — 2ie 
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and the third to 



2i5rs I dydzdudv 



duo I duj2 

l/^G,,(y,u)(-2V:r)G.,(u,v)(-^V:r)Ga.,(v,z)/G^,(z,y) 

1 1 



(4.13) 



-zB f ,t-B 

dt—^ 

B 



t — ul)i+ ul)2 — 2ie t — UJ2 + uji — lie 



The next step is to consider separately the expressions arising from the two terms 
in the last factor. For the first term, close the uj\ contour up and for the second 
close the UJ2 contour up. Here the asymmetry between uj\ and 0^2 causes some 
differences from the corresponding simplification of the first diagram. The first 
terms involve both a double pole and and three simple poles at uj\ = —B + ie. 
Whereas the second terms involve only a single simple pole. The four simple pole 
contributions from each diagram lead to expressions very like (4.11). Namely, 



t-B. 



-B5rs I dt—^[{B\y'' ^ -"^ 



^ p-G_,_5p--J-^^^- \B) 



B 



p' + B 



p' + B' 



^ p'"G -By' -yr^y' \B} 



<^l'''"? + B + < 



p' + B' 



+ {B\p'- 



y^G^t-BV^ 



p^ + B + t 



p^ + B + t 



p"^ \B) 



\B\p-^^^y^\B){B\y'' 



k 1 



p'^\B) 



+ (^1 y'zyj^y^'G-Bf^ ^,, , ^„ , v"^ m 



p2 + 5^ 



p^ + B + t' 



(4.14) 
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for the first diagram, and 



oo 



B 



1 1 

p^ + l3 + t p^ + l3 



+ ^B\p^ 2 , o , . f^G^BV^'^rv^y^ \B) 



(4.15) 



- 7 {B\ y^^^p"" \B) {Blp'^^^y' m 

t p^ + i3 p^ + 13 



for the second. (4.14) and (4.15) do not include the double pole contribution to the 
integral of the first terms over uji. These arise from the ground state contributions 
to the second and fourth G"s in (4.12) and to the first and third G"s in (4.13). 
Evaluating the residue of the double pole in uji and then doing the uj2 integral 
leads to the double pole contributions 



oo 



„, , , t- B d 

BOrs I dt 



t2 du 
B 

oo 



^^^'«'^^i^j^^^ <^i'''" p^+fiV.-y "i^> 



Jcj=0 



B 

(4.16) 
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and 



BSr'i / dt 



t-B d 



t2 du 



{B\y' 



B 



p^ + B -CO 



/J- ID 
dt ^ 

B 



[{B\y' 



(p2 + Bf 



y'\B){B\p'' 



y'\B){BW 



B + t-u 



p'^\B) 



u}=Q 



p2 + 5 



y'\B){B\f 



1 



(p2 + 5 + t)2 



p2 + 5 + t 

p- 15)] 



p-|5) 



(4.17) 



respectively. 

Notice that lines 2 and 5 of (4.14) and lines 4 and 6 of (4.15) cancel lines 2, 4, 
6, and 8 of (4.11) by virtue of the fact that 

{B\ y^^^p""' \B) = - {B\p-^^]—y^ \B) 

p^ + ±5 p^ + ±5 

by time reversal invariance. This cancellation together with similarity of structure 
make it helpful to regroup the various contributions by gathering together all terms 
with the central Green function G^t-B in group 1, 



-B5r. I dt-^\{B\y^ 



t2 



B 



+ {B 
+ {B 
+ {B 
+ {B 
+ {B 
+ {B 
+ {B 



p2 + fi 



P' 



p'^G^t-By' 



1 



p2 + 5 + t 



p-^\B) 



^ y^G_,_^p-^-l-^y^ \B) 



p^ + B + t 
1 



P __2 ■ n ■ - ^ G^t-BP 



p2 + 5 

. 1 



y 



p^ + B + r 
■^^^P^G^t^By''- 



p2 + fi 

1 



y' \B) 



p2 + fi^ 



2 + fi + t 



P 



p- 15)] 



y'l^P"'G^,^BP'^^^y^\B) 

p"" , , \ , y G^t-BV 



p' + B' 



k 



-p"^ \B) 



p' 






y'TI^P'^G^t^BP'^Zy^y' m, 



p2 + 5^ 



p2 + 5- 



(4.18) 
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and all those with central Green function G^b into group 2, 



oo 



B 



kr" ^ 



Y> + B + t p^ + -D 

The double pole contributions are also naturally associated with group 2. 

In the appendix we describe the detailed evaluation of all these expressions. 
By far the largest contribution comes from the second and third diagrams of 
Fig. 3. The total contribution of the first diagram is approximately — 0.055197/i?^ 
or about 16.6% of the random phase approximation. Compare this to the total 
contribution of the second and third diagrams of +0.50913/5^ or —153% of the 
RPA. In the next section we show that this huge correction can be attributed to 
a substantial shift in the ionization energy of the open polymer. The large terms 
in perturbation theory can be identified as "self energy" corrections which can be 
easily summed as a geometric series in the usual way. 

5. Summing Large Corrections 

We have found very large corrections to RPA in next to lowest order of per- 
turbation theory. But we have also seen that all the large contributions come from 
diagrams that have the interpretation as self energy bubbles on the propagators 
entering the RPA. Thus we can hope that an improved RPA would result if we 
use the full propagators rather than the bare ones. In other words, we should take 
the sum of diagrams shown in Fig. 4 as the propagator to use in RPA. An even 
more sophisticated improvement would be to replace each bubble in Fig. 4 with 
the RPA sum of Fig. 2. 
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Fig. 4. Summing self energy diagrams. 

More generally we can define a "self energy" part, S, so that iS is the sum of all 
the (amputated) one particle irreducible diagrams contributing to the correlation 
function of ip and ip'. Then the full propagator 

{G\T[My,t)4{^,0)]\G) ^ -t6r,sA 

satisfies an integral equation 

A = G + GEA. (5.1) 

Because the operator ipr annihilates a bond in the polymer, the location in —uj of 
the lowest pole in A gives the "exact" ionization energy, which in zeroth order of 
perturbation theory is just B, the lowest pole in G. To find this lowest pole it is 
sufficient to look at a matrix element of A between any two states which couple to 
it. Assuming that the zeroth order ground state continues to couple to this pole, it 
is convenient to use a new concept of "irreducible" to mean that one can't cut the 
diagram in two by cutting a single line in its ground state. This new irreducible 
part will then be called S'. A formula for S' is obtained by splitting G = G + Gq 
(see (4.8)) where, in frequency space, 

^ _ \B){B\ 



—B — u + ie 
Then we have the integral equation 



E' = S + SGS'. (5.2) 
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That is, we put back into S' all those reducible contributions for which the single 
line is in an excited state. 

With these definitions a little algebraic gymnastics leads to the following for- 
mula for the full propagator: 

A = G + Gi:'G + iI + Gi:')—- ^—-j -PoiI + i:'G) (5.3) 

—B — uj — Wq + te 

where Pq = \B) {B\ is the projector onto the zeroth order ground state, and Sq = 
{B\Ti'\B) a real number. As we explain in the Appendix, the first step in the 
evaluation of a typical diagram is to close one or more of the frequency integrals 
in the upper half plane, picking up some of the ground state poles in the various 
G^j's. In a more sophisticated application of perturbation theory, where A is used 
in place of G, this step will pick up the pole at u = ujq where 

D{uJo) = -B-uJo-i:'Q{uJo)+ie = 0. (5.4) 

The residue of this pole is 

^ {I + G^,^'M)Po{I + ^'MG^,) = —L_G^^^'{ujo)PoJ:'{ujo)G^,. (5.5) 



Note that the form on the r.h.s. would be meaningless if ujq = B since G-s is 
infinite. But since the shift of ujq away from its zeroth order value is taken into 
account, this is not a problem. (If by some fiuke the shift were zero, then S'_^ 
would supply a zero multiplying the infinity, leaving the r.h.s. indeterminate.) 

Now we can use these results to specify our improved version of RPA. It is 
obtained by replacing each of the bare propagators in the first diagram of Fig. 2 
with the the full propagator (5.3). Corrections which involve "self energy" inser- 
tions such as the second and third diagrams of Fig. 3, will then automatically be 
included in this approximation. The first diagram of Fig. 3 is a typical diagram 
not included. Our experience with straight perturbation theory indicates that the 
diagrams not included give roughly a 15% correction. 
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Inserting (5.3) into the RPA diagram, we can carry out the uj integral, just as 
with the bare propagators, by closing the contour into the upper half plane, picking 
up the lowest pole at cj = ujq, determined by (5.4). The +ie of (5.4) amounts to 
the contour prescription which states that this lowest pole is the only singularity 
picked up by this procedure. The result is 

TTimproved,jfc/n^i _ _9A (-5| T^' {ujQ)Gu,oy^[G^,, + G ^,,1:' {ujq)G ^,]y G^JI'Jujq) \B) 

(5.6) 
Note that the pole contribution to A vanishes when sandwiched between y^ and y 
by rotational invariance. If we were given the exact result for S' we could use this 
formula to calculate the string tension. We would first have to find cjq by solving 
(5.4). Then evaluating S' at u;o and plugging into (5.6) would give us the required 
information. 

Unfortunately, S' is not a quantity we can evaluate exactly, and to assess the 
ramifications of this procedure we must make further approximations. A plausible 
approximation would be to calculate S in perturbation theory. Thus in lowest 
order we would simply take the second diagram of Fig. 4, with external propagators 
removed, for S(c<j). One could then try to use this approximate form for S in (5.2), 
to determine S'. Although inexact, this is a substantial improvement over straight 
perturbation theory because, for instance, the procedure does not violate easily 
proved positivity constraints on Qrs- Probably the approximation most consistent 
with RPA would be to use the RPA sum for the bubble in the lowest order diagram 
for S, as we described in the opening paragraph of this section. In the remainder 
of this section we attempt a simplified version of the first mentioned scheme. The 
complete implementations of the first scheme and of the second would involve some 
nontrivial technical developments, and we leave them for future work. 

To get a feeling for how these considerations affect the numerical results of 
the previous section, we close this section with the evaluation of the string tension 
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using the lowest order approximation to S, 

S(cj) ^ 

B f q^ I sr^ p^\m){m\p^ p^|0)(0|p^ 

aq- 



TT y ^(q2 + 5)2 I ^ E^ + q2 + 5 _ ^ _ 3ie Eq + Ci^ + B - uj + ?>ie 
\m>0 / 

in an approximate version of (5.2). Recalling that the large corrections found in 
the previous section could be traced to the contribution of the second term of (4.7), 
it would seem to be reasonable to truncate the G appearing in (5.2) to that term. 
We therefore define, in notation explained in the appendix, 

^ - '*'™«(^' - BH-.IB) !--> <--! - ^:hB i«> <«i ■ '"' 

We can then use (5.2) to determine by straightforward algebra the matrix elements 

(i?| S' |i?), (— cj| S' |i?), {B\ S' |— cj), and (— cj| S' \—uj) in terms of analogous matrix 
elements of S. 

Define 7^ = l/{B\n{-^ / B)}, 70 = l/(cu + 5), and 

^1 = 1-7.; (-^1 S l-cc.) - 70 (5| S |5) + 

7^7o[(5| S |5) (-cu| S l-cu) - (-cu| S \B) {B\ S |-a;)]. 



Then we find 

{B\^'\B) 
{-u\T.'\B) 

{B\T.'\-uj) 
(-cj|S'|-cj) 



[B\ S |5) - 7^[(fi| S |fi) (-cj| S |-cj) - {-u\ S |fi) (fi| S \-u)] 
(-cj|S|5) 






-uj\ S |-cj) - 7o[(B| S |B) (-cj| S \-uo) - {-uo\ S |5) {B\ S |-cj)] 



For the matrix elements of S, we find, after doing angular integrals and some minor 
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changes of integration variables, 



oo oo 



(B\^\B)=Bi^[^ ^^::'^<r'L ^BfM 



v^ J u'^ u + V — I — uj/B 
1 1 



(5|S|-cj) =B 



oo oo 

dv f du 



{v-l){u-l) 



v^ J u {u-1- uj/B){u + V-1- uj/B) 
1 1 



= 5/2(p) 



-uj\^\B)={B\^\-u) = Bf2ip) 

du 



(-cu|S|-cu) =B 



oo oo 

dv 



{v-l)iu-l) 



v"^ J {u-l-uj/BYu + v-l-oj/B 
1 1 



= Bh{p) 



where we have defined p = — 1 — uo/ B. The v integral is elementary, with the 
results 



(u — l)du 



flip) 



/2(P) 



Mp) = (u- l)du 
1 



ln(l + M + p)-l ln(l+M + p)" 




ln(l+M + p)-l ln(l + M + p) 



{u + pY 



{n + pY 



Putting everything together, we express all quantities in terms of the functions 
/i, /2, h- For Di we have 



D^ ^1 



h{p) , Hp) h{p)h{p)-f2{pf 



ln(l + p) 



P 
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pln(l + p) 



and for the matrix elements of S', 



{B S' B) ?: 




{-u S' B) P 


. r/2(p) 


{B\^'\-uj) ^ 


. r/2(p) 


(-u\J:'\-u) p 


r, f sip) + [flip) f sip) -f2ip?]/p 



Using the above approximate resuh for {B\ll' \B) = Sq in (5.4) gives the following 
approximate equation for po = — 1 — ^o/B: 

Mpo) - [fliPo)hiPo) - /2(po)Vln(l + po) ^ 
^l(Po) 

Using a simple iterative search for solutions to this equation leads to a value for po 
between 0.589941 and 0.589942, so we can say that po ~ 0.590. This corresponds 
to a shift in the ionization energy from B to 1.5905. The same numerical data 
give 1 + dTiQ/duj ^ 1.06, which is consistent with unity in the present discussion. 
Straight perturbation theory would lead to the result uq = —(1 + /i(0))i?. But 
/i(0) = 111 — /21 ~ 0.59086. So, curiously enough, the improved prediction for uq 
is barely distinguishable from the calculation of the first two orders of perturbation 
theory! As an aside, we recall that in the Hartree-Fock method which we are using 
here, the relationship of the shift in ionization energy (or single particle energy) is 
related to the shift in the ground state energy per bond (or per single particle) by 
a factor of 2. Thus this 60% shift in the ionization energy, which is rather huge, 
corresponds to a more moderate 30% shift in the ground state energy. 

Having approximately determined ujq, we consider the evaluation of (5.6) at 
this value. In the same spirit as our approximate treatment of (5.2), It is reasonable 
to approximate the first and last G^o's on the r.h.s. by 

^-"i?ln(-c.o/5)'"''°^^"''°'- 
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Then 

■rTlmproved,jfc/r|\ ^ 
yy,rs v"/ ^^ 

Referring to our experience with straight perturbation theory the second term 
includes contributions in next to lowest order which are relatively small, i.e. roughly 
15%. Thus it is fair to drop it in the present approximation. The resulting matrix 
element is easily evaluated 



652(1 + po)3 ■ 6fi2- 
Thus to an expected 15% accuracy we have 

nlmproved,ifc(o) ^ -'^^rs'—^f^^ {-U^olv' G^oV' h^o) ^ -0.333^. 



Thus our improved RPA causes predicts a reduction in 11 by approximately a factor 
of 3 over the naive RPA. But now we have explicitly incorporated all of the large 
corrections in next order in a partial resummation of perturbation theory. The 
effects left out amount to a roughly 15% correction in this next order, and it is 
reasonable to expect that this is a rough indication of the accuracy of our improved 
RPA. The factor of 3 reduction in 11 corresponds to a factor of \/3 increase in the 
rest tension leading to 



nlmproved RPA 



1 



~ ^./Improved RPA ~ "^ 
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6. Conclusions 

In this article we have initiated a detailed quantitative study of the "wee par- 
ton" formulation of string theory proposed in Ref.[3,9]. Although the quantitative 
relationship between the microscopic string bit dynamics and the low energy stringy 
properties of string theory is subsumed in merely a renormalization of the string 
tension, we believe that these details are extremely important for understanding 
any processes which are capable of probing the microscopic structure of the theory. 
A gedanken experiment which could in principle reveal this structure might be the 
physics near the horizon of a black hole as perceived by an observer far away. 
It is, of course, unlikely that the specific string bit model we have chosen to analyze 
in this article is the ultimate underlying dynamics of the bosonic string. Our pur- 
pose, rather, is to pick a model and explore its detailed dynamics using well-known 
methods in many body physics. The model we have chosen, with its delta func- 
tion interaction potential, does nonetheless have several appealing features which 
recommend it, including locality in the transverse space and (anomalous) scale 
invariance. 

For pedagogical simplicity we have only tried to account for the bosonic string 
in this article. But one reason for subjecting the string bit dynamics for this 
unrealistic model to such scrutiny is that the superstring could be much more 
sensitive to the microscopic details than is the bosonic string. In that case we need 
to have a high degree of confidence in the techniques used to analyze a superstring 
bit model. We hope to devise such a model in the near future. 

The numerical results of our study are encouraging. Although the most naive 
application of RPA to the calculation of the rest tension is subjected to huge 
(150%) corrections in next order of perturbation theory, we have shown that the 
bulk of these corrections can be incorporated in a Hartree-Fock like improvement 
of the zeroth order RPA. The corrections that remain are estimated to be around 
15%. We did make justifiable simplifications of the Hartree-Fock treatment, but 
there is certainly room for improvement. It would be particularly desirable to go 
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beyond the lowest order approximation for S. Also, in calculating S' we truncated 
the Green function G' by dropping the "free" contribution which was shown to be 
relatively small in perturbation theory. However, it should be a tractable problem 
without this truncation. Thus, although the simplifications we made were consis- 
tent with our approximation, we expect that they can be avoided with a little more 
effort. 
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APPENDIX 

In this appendix, we describe the evaluation of the various matrix elements and 
integrals occurring in (4.18) and (4.19). The Green function occurring in group 1 
is not singular in the integration range, so all e's may be set to zero. Then there 
is no distinction between G and G, so referring to (4.7) and going to momentum 
representation, we may write 

r 1 I \t + B){t + B\ 

where we have introduced the state |t + B) defined by 



B 1 



TT p^ + i? + t 

Defining 
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we can also write 



G. 



1 



B 



P' 



B 2B 



\B){B\ + \B'){B\ + \B){B'\ 



(A.2) 



We first consider the contributions due to the first terms in (A.l) and (A.2), 
the "free" Green functions. The first four matrix elements in (4.18) turn out to 
equal each other. Taking the first one, for example, we get 



1 



d 



p" 



B f 2ip^p^ 

77 P (p2 + Bf (p2 + 5 + tf'^ (p2 + 5 + t)(p2 + B) 



5k jB / du 



B 



252 AB 2 



u^ 



u^ V? {u + t)3 



252 3B 1 ■ 

u^ -u^ u^ 



[u + ty^ 



Let us define 



v-1 



J-n,m — / f^"^' 2 



du 



1 



u^{u + vy 



(A.3) 



Needed numerical values for some of the In,m are displayed in the following table. 



^n,m 


n\ m 


1 


2 


3 


4 


5 











0.007530 


0.001851 


1 


0.886294 


0.090863 


0.017005 


0.004129 


0.001163 


2 


0.295431 


0.045431 


0.010202 


0.002751 


0.000832 


3 


0.170431 


0.029519 


0.007168 


0.002037 




4 


0.118925 


0.021702 


0.005493 






5 


0.091148 


0.017114 








6 


0.073836 
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Then the "free" parts of the first four hnes of (4.18) evaluate to 



^"^-[2/4,3 - 4/3,3 + 2/2,3 + 2/5,2 - 3/4,2 + h,2] ~ ^o ' • (A.4) 



The 5th and 8th matrix elements are also equal, and their free contributions to 
(4.18) are easily shown to give 



-^SrsSjk.. ^j , . -(U)41865^ 

— g2 — [J6,i - 2^5,1 + M,iJ ~ ^2 • y^-^) 



The 6th and 7th matrix elements in group 1 are the most tedious to evaluate, but 
the steps are straightforward, with the result for their free contribution to (4.18) 



or r 

^f^ [2/2,5 - 4/1,5 + 2/0,5 + 4/3,4 - 6/2,4 + 2/1,4 + 2/4,3-2/3,3 + /2,3] ~ 

-0.0149325„5jfc 

(A.6) 
Finally the free contribution to (4.19) is the same for each of the 4 matrix elements, 
leading to a total contribution 



■[-6/6,1 + 10/5,1 - 4/4,1] ~ So -■ (A-7) 



We record here the total contribution from the "free" term of the Green function 
to (4.18) and (4.19), the sum of (A.4), (A.5), (A.6), and (A.7), -0.033284/52, 
which is approximately 10% of the lowest order (RPA) term. 

Next we turn to the contribution from the second term of (A.l) to the various 
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terms in (4.18). The evaluation requires the matrix elements 



{t + BW , \,. P'^\B) =^^^ [dp- P' 



p2 + 5 + t^ I / ^ J ^(p2 + 5)(p2 + 5 + ^)3 



oo 



^6,m f^^ u-l 



B J u{u + t/Bf' 
1 

{t + 51 P^^^y^ \B) =ZlEhiL f dp- P' 



p2 + fi^ ' ' TT J ^(p2 + fi)3(p2 + 5 + t) 



CX) 



-.5,^ /^^ .-1 



E J u^{u + t/By 
1 



and their complex conjugates. The integrals here are elementary: 



oo 

, M-1 1 /111 1,, , 

du— -^ =— — - + - - -(1 + -) ln(l + v] 

1 

, «-l 1/11 I,, 1,, ,, 

du^^-. =-- + -- - 1 + - In 1 + V, 

u-^yu + v) V \v I V V 

1 



Using these results we find that the second term of (A.l) inserted into (4.18) yields 

45^ r t-l /I 1 1 1 y ^^ 0.005795^,,^ 

(A.8) 
for the first 4 lines of (4.18). For the 5th and 8th lines of (4.18) the corresponding 
contribution is 



00 „ 

-2^r.^,fc [ t-l (I 1 1 1 y ^^ -0.009937^,,^,fc 

52 y t4in(l + t)Vt 2 t^^+t^^^^^ + ^^j -- 52 

(A.9) 
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and for the 6th and 7th hnes, 



-25rR5 



rs"jk 



dt- 



t-l 



1 



-0m09676r.S 



rs"jk 



52 

(A.IO) 
The numerical values for these integrals were obtained by a simple integration 
routine. For reference purposes we also quote the values of the terms in (4.11) 
that have been cancelled, so that the total value of the first diagram in Fig. 3 can 
be worked out. This works out to a total of —5rs6kj/{^^B'^). Noting that (A. 8) 
almost cancels (A. 4), we see that this last bit is the dominant contribution to the 
diagram, about 16% of the RPA contribution. 

The contribution from the last three terms of (A. 2) to (4.19) requires the matrix 
elements 



{B\y 



y \B) — 



p 



2 + 5- 



(51/ 



1 



6fi2 

6 



5fi2 



' ^in 



[B' p 



'— p^ \B) = I du^ — \— 

1 

00 

—^ p- \B)=-- I du ,,''~\ , . 

p2 + 5 + r ' ^ Bj u^{u + t/B) 

1 



Using these results, we then obtain for the remainder of the evaluation of (4.19) 



-4:6rsS 



rsOjk 



52 



— - ^)(-^l,l ^ -^2,1) - g(-^2,l - -^3,1) 



5r.<i5 



rsOjk 



■7 2 

Y^(/i,i - h,i) + 2(/2,i - /3,l) 



0.359069 
fi2 



(A.ll) 
Finally, we evaluate the double pole contributions (4.16) and (4.16), which requires 
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the additional matrix elements 



1 fc I D\ ^Jk 



{B\y \ , ^^y \B) 



oo 



p2 + 5 + t)2^ ' ' 5 7 u\u + t/B)^' 

1 



Then the total double pole contribution evaluates to 



52 



-(/l,l-^2,l) + -(/l,2-/2,2) 



0.133316 
' 52 



(A.ll) 



The grand total of all contributions is 

TOTAL 



A539925r s5kj 
52 ^ 



which is —136% of the RPA result. Clearly this cannot be regarded as a small 
correction! 
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